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Abstract 

We study the absorption of massive spin-half particles by a small 
Schwarzschild black hole by numerically solving the single-particle Dirac 
equation in Painleve-Gullstrand coordinates. We calculate the absorption 
cross section a(E) for a range of gravitational couplings Mm/mp 2 and 
incident particle energies E. At high couplings, where the Schwarzschild 
radius Rs is much greater than the wavelength A, we find that a(E) 
approaches the classical result for a point particle. At intermediate cou¬ 
plings, Rs ~ A, we find oscillations around the classical limit whose precise 
form depends on the particle mass. These oscillations give quantum vi¬ 
olations of the equivalence principle. At high energies the cross section 
converges on the geometric-optics value of 277rf?|/4, and at low energies 
we find agreement with an approximation derived by Unruh. When the 
hole is much smaller than the particle wavelength we confirm that the 
minimum possible cross section approaches 7rJ?|/2. 

PACS numbers: 04.70.Bw, 03.80.+r, 04.62,+v, 03.65.Pm 


1 Introduction 

It is widely accepted that general relativity and quantum mechanics are incom¬ 
patible in their current form, yet a theory reconciling the two remains elusive. 
Despite this problem, it is only at the smallest length scales (l < Ip), or highest 
energies, that we expect substantial modification to existing theory. At low 
energies, away from spacetime singularities, the propagation of quantum fields 
on gravitational backgrounds is well understood (see the books by Birrell & 
Davies JT| or Chandrasekhar 2 ]). 

Interest in the absorption of quantum waves by black holes was reignited in 
the 1970s, following Hawking’s discovery that black holes can emit, as well as 
scatter and absorb, radiation [3|. Hawking showed that the evaporation rate 
is proportional to the total absorption cross section. More recently, absorption 
cross sections (or “grey body factors”) have been of interest in the context of 
higher-dimensional string theories. 

In a series of papers 011 El Sanchez considered the scattering and ab¬ 
sorption of massless scalar particles by an uncharged, spherically-symmetric 
(Schwarzschild) black hole. Using numerical methods she showed that the total 
absorption cross section (as a function of incident frequency) exhibits oscilla¬ 
tions around the geometric-optics limit characteristic of diffraction patterns. 
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Unruh '7J studied the absorption of massive spin-half particles by piecing to¬ 
gether analytic solutions to the Dirac equation across three regions. He showed 
that, in the low-energy limit, the scattering cross section for the fermion is ex¬ 
actly 1/8 of that for the scalar particle. He also derived an approximation to 
the total cross section valid at low energies, which we revisit in sectional 

In this paper we return to the massive fermion absorption problem studied by 
Unruh. We employ a different coordinate system, but retain equivalent ingoing 
boundary conditions at the horizon. Instead of using analytical approximations 
we numerically integrate the Dirac equation to calculate the absorption cross 
section over a range of energies and gravitational couplings. We compare our 
results with the classical cross section for a massive particle, and with Unruh’s 
low-energy approximation. 

The natural dimensionless parameter to describe the strength of the gravi¬ 
tational coupling between a black hole (of mass M ) and a quantum particle (of 
mass m) is given by 

GMm Mm 7 tRs . . 

a = —r- = — 2 “ =-r— (1) 

nc m P Ac 

where Rs is the Schwarzscliild radius of the hole, Ac is the Compton wavelength 
of the quantum particle, and m p is the Planck mass. We use the symbol a 
because it has an analogous role in gravitation to the fine-structure constant 
in electromagnetism. We expect quantum effects to be important when a ~ 1, 
whereas in the high-a limit classical effects should dominate. 

In first-quantised theory the capture of light and matter by a black hole 
is a one-way process. The direction of time implied by this process is not 
revealed in Schwarzschild coordinates, however, as these are manifestly time- 
reverse symmetric and are invalid at the horizon. Time-asymmetric coordinates, 
such as Eddington-Finkelstein coordinates, allow the continuation of the metric 
across the horizon and allow us to correctly study the properties of wavefunc- 
tions El El EDI- We then find that ingoing states correspond precisely to those 
that are regular at the horizon. 

Here, rather than using Eddington-Finkelstein coordinates, we prefer to 
work with the coordinates first introduced by Painleve EH and Gullstrand Ga¬ 
in these coordinates the metric becomes 

ds 2 = f 1- j dt~ — J - dt dr — dr 2 — r 2 dfl 2 (2) 

The utility of this form of the Schwarzschild solution has recently been high¬ 
lighted by Martel & Poisson m and others 0- For black holes (as opposed 
to white holes) the negative sign for the crossterm dtdr is the correct choice, 
as this guarantees that all particles fall in across the horizon in a finite proper 
time. This sign is also uniquely picked out by models in which the black hole 
is formed by a collapse process {HJ. One advantage of this system is that the 
time coordinate has a natural interpretation as the proper time measured by an 
observer in freefall starting from rest at infinity. 

According to general relativity the classical absorption cross section of a 
Schwarzschild black hole is given by 

aabs = Ij^T ( 8m4 + 20u2 “ 1 + (! + 8lf2 ) 3/2 ) ( 3 ) 
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where u is the velocity of the particle j7j. In accordance with the equivalence 
principle, the cross section is independent of the particle mass m. We expect 
that the quantum cross section will approach this value in the limit a>l (that 
is, Rs > A). 

We start with the radial separation of the Dirac equation in Painleve- 
Gullstrand coordinates. We then study the properties of solutions around the 
horizon, identifying the physical branch of regular solutiuons. For unbound 
states E > me 2 we find that the physical solutions are composed of ingoing and 
outgoing waves at infinity. By numerically finding the ratios of these waves in 
any given angular mode we are able to compute the absorption spectrum. We 
use natural coordinates G = Ti = c = 1, except in cases where inclusion of the 
factors adds clarity. 


2 The Dirac equation 

We let { 70 , 71 , 72 , 73 } denote the gamma matrices in the Dirac-Pauli represen¬ 
tation, and introduce spherical polar coordinates (r, 9, <j>). From these we define 
the unit polar matrices 

7 r = sin$(cos<ji>7i + sin</>72) + cos$73 
70 = cos0(cos<(>7i + sinc^72) — sin$73 
'fy = — sin 0 71 + cos (j> 72. ( 4 ) 

In terms of these we define four position-dependent matrices {fft, g r , ge, g<t>} by 

[2~M 

9t= 70 + y —' 7r 
9r = lr 

These satisfy the anti-commutation relations 

{9(1,9*} = 2 g^ 1 

where g is the Painleve-Gullstrand metric of equation (J5J). The reciprocal 
matrices {g*, g r , g e , g^} are defined by the equation 

{<A<?„} = 2^J, (6) 

and both sets are well-defined everywhere except at the origin. 

The Dirac equation for a spin-half particle of mass m is 


ge = r'yg 

g<t> = r sin 6 » 70 . (5) 


ig^Vpip = mif), (7) 

where 

£a/9 = £[7a,7/j]- ( 8 ) 

The components of the spin connection are found in the standard way Ha and 
are particularly simple in the Painleve-Gullstrand gauge (10) . 


■"-r^Sa/s =- -i- \l — 


4 r 


2 M 


-7o- 


( 9 ) 
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An advantage of our choice of metric is that the Dirac equation can now be 
written in a manifestly Hamiltonian form 


iftif) — *7° 





mi/j, 


( 10 ) 


where ft is the Dirac operator in flat Minkowski spacetime. The interaction 
term is non-Hermitian, as the singularity acts as a sink for probability density, 
making absorption possible. 

The Dirac equation is clearly seperable in time, so has solutions that go as 
exp (—iEt). The energy E conjugate to time-translation is independent of the 
chosen coordinate system, and has a physical definition in terms of the Killing 
time D2H- We can further exploit the spherical symmetry to seperate the spinor 
into 




-™( urWxKCM) A 

r Vu 2 (r)ovx£(M)/ 


( 11 ) 


where 


oy = sin(9(cos()><7i + sin^>cr 2 ) + cos 0 tJ 3 . 


( 12 ) 


The angular eigenmodes are labeled by k, which is a positive or negative nonzero 
integer, and /i, which is the total angular momentum in the 9 = 0 direction. 
Our convention for these eigenmodes is that 


(a-L + h)x£ = Khx£, K = ±(j + 1/2) = ...,-2,-1,1,2,.... (13) 


The positive and negative n modes are related by 


(JrX^ = X- K 


(14) 


and are normalised so that 

J d(j> J dd sin 9 x£(6>,<(>) t x£,(0,0) = L'V' ( i5 ) 

The trial function GB results in a pair of coupled first-order equations 


(2) "Ur T) 

n/r i{E + m) — (2M/r) 1 ^ 2 (4r)~ 1 

i(E — m) — (2M/r) 1 / 2 (4r) _1 — n/r 

The equations have regular singular points at the origin and horizon, as well 
as an irregular singular point at r = oo. As far as we are aware, the special 
function theory required to deal with such equations has not been developed. 
Instead we use series solutions around the singular points as initial data for a 
numerical integration scheme. 



3 Series Solutions and Boundary Conditions 


As is clear from lP|) . there is a regular singular point in the coupled equations 
at the horizon, r = 2 M. We look for series solutions 


U = 



(r 



(r - 2 M) k 


(17) 
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where s is the lowest power in the series, and a*,, bk as coefficients to be deter¬ 
mined. On substituting into m and setting r = 2 M we obtain an eigenvalue 
equation for s, which has solutions 

s = 0 and s = — \ + 4 iME. (18) 


The regular root s = 0 ensures that we can construct solutions which are finite 
and continuous at the horizon. We will see later that regular solutions auto¬ 
matically have an ingoing current at the horizon. The singular branch gives 
rise to discontinuous, unnormalisable solutions with an outgoing current at the 
horizon [5]. We therefore restrict attention to the regular, physically-admissable 
solutions. The eigenvector for the regular solution has 

( a o\ _ (K-2iM(E + m) + l/4\ , , 

\b 0 J \K + 2iM(E-m)-l/iJ ' [ ’ 

In order to expand about infinity we need to take care of the irregular 
singularity present there. There are two sets of solutions, t/*° ut ) and U^ ln \ 
which asymptotically resemble outgoing and ingoing radial waves with addi¬ 
tional radially-dependent phase factors. To lowest order, 


[/(out) 


[/M 


_ ipr i(4>i+4>2) | 1 I 

{p/(E + m)J 

— p -ipr p i{4> l-0 2 ) | 1 | 

{-p/(E + m)J 


( 20 ) 


where the phase factors </>i (r), cj> 2 {r) are given by 

0 i(r) = EV8Mr, fair) = — (to 2 + 2 p 2 ) ln(pr) (21) 

p 

and the momentum p is defined in the usual way, p 2 = E 2 — to 2 . The general 
(regular) solution as r —» oo is a superposition of the ingoing and outgoing 
waves, 

U(r oo) = a K U^ + /3 K I/ (out) (22) 

for each angular mode. The magnitudes of a K and f3 K determine the amount of 
scattered and absorbed radiation present. 


4 Absorption 

The spatial probability current is conserved for states with real energy, E > to. 
For each angular eigenmode we obtain a conserved Wronskian W K 

W K = (u\u\ + u\u 2 ) — y/2M/r(uiu\ + U 2 W 2 ) (23) 

which measures the total outward flux over a surface of radius r. At the horizon 

W K = -\u 1 -u 2 \ 2 (x-\a 0 -b 0 \ 2 (24) 

so the flux is inwards for all regular solutions. On substituting the asymptotic 
forms of equations CT and lt22ll into equation l|23|) we find an expression for 
the Wronskian in the large-?’ limit, 

Wn = -=^-(\a K \ 2 ~m 2 ) (25) 

E + to. V / 
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The coefficients a K and /3 K can be determined (up to an overall magnitude 
and phase) by matching the ingoing (regular) solution at the horizon to the 
asymptotic form in the large-r limit. We choose the normalisation of each 
angular mode so that W K = —1, and write the most general solution to the 
wave equation as 

->P = 'Y1 i 9 k ,^k (26) 

Mo 

where ijj K are spinors of the trial form with ui, 112 as given by CT . 42U, and g K 
are complex coefficients. The total absorbed flux is then just 

w t ot = Y,\g« l 2 - ( 2 ?) 


We now employ a partial wave analysis to derive a simple expression for the 
absorbed cross section. We write the asymptotic behaviour of il> as the sum of 
a plane wave (propagating in the 9 = 0 direction) and an outgoing scattered 
wave, 

^ = e i P rc OS e^ i+ /W e i P r^ 9 (28) 

r 

where d' 1 , 4'2 are constant spinors. The plane wave can be decomposed into 
ingoing and outgoing radial waves in the large-r limit. We equate the ingo¬ 
ing part of the plane wave with the ingoing part of the asymptotic wave (l22ll . 
Normalising the plane wave to 2 E particles per unit volume we find 


= i(-l r+ WME + m) , 

2 p /N 

for each angular mode. The total absorption cross section <r a b s is the ratio of 
the ingoing flux am to the incident flux of the plane wave (2 p), 


_ 71 N 

abs 2p{E-m) ^ \a K \ 2 ' 


(30) 


At low energies, the ft = 1 states dominate the absorption, but at higher 
energies we need to sum over a range of ft. 


5 Results 

We determine the coefficients a K required for m by matching the ingoing 
solution at the horizon to the asymptotic form |22|) at infinity. In a similar 
calculation, Unruh |7J used analytic approximations to the radial functions to 
find the leading contributions to the cross section. Here, we use numerical 
integration of the wavefunction out from the horizon to match the solutions, 
and compare our results with the analytic approach. 

Figure 0 compares the result of our matching calculation at a = 0.2 with 
the classical cross section of a point particle © over a range of energies. We see 
that the quantum absorption cross section oscillates around the classical value, 
as found by Sanchez 00 for the massless scalar wave. For a given a we find 
the period of these oscillations is approximately constant, but the amplitude 
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Energy, E / me 2 

Figure 1: Classical and quantum absorption cross section. The plot com¬ 
pares the absorption cross section for the Dirac wave [solid] with the clas¬ 
sical prediction for a point particle [dotted], for a gravitational coupling of 
a = Alm/nip 2 = 0.2. The cross section is plotted in units of ( GM/c 2 ) 2 (pro¬ 
portional to the event horizon area), and the energy in units of the rest mass 
energy me 2 . 
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Figure 2: Quantum absorption cross section for a range of couplings. The plot 
shows the cross section at a range of gravitational couplings, 0.05 < Mm/m p 2 < 
0.2. The horizontal line is the photon limit. 


decays as E —> oo. As a is increased we find that the magnitude and period 
of the oscillations decreases. In the a 1 limit we recover the classical cross 
section. 

Figure Q illustrates that the precise form of the oscillation depends on a, 
and therefore on the mass of the quantum particle. This represents a quantum- 
mechanical violation of the equivalence principle. At sufficiently high energies, 
we see that all cross sections tend to the photon limit of cr a b s = 277 t{GMc~ 2 ) 2 . 
As noted by Unruh, all particles travelling close enough to light speed, u ~ 1, 
see a black hole of roughly the same size, regardless of mass or spin. 

Unruh also showed that in the low-energy limit, the Dirac cross section is 1/8 
of the scalar cross section, and absorption is dominated by the lowest angular 
momentum modes, |k| = 1. In this limit he showed 

Jabs _ _47 t 2 (1 +u 2 )a _ 

[GMc~ 2 ) 2 u 2 ( 1 — m 2 ) 1 / 2 {l — exp(—27ra(l + u 2 )/u{ 1 — u 2 ) 1 / 2 )} 


This proves to be an excellent fit to the numerical cross sections in the low-energy 
regime, such as those shown in figure 0 A minimum-possible cross section can 
be found by considering the a ~ 0 limit of equation which reduces to 


^abs 27T 

( GMc - 2 ) 2 ~ ~u' 


(32) 


The minimum value of this occurs in the u = 1 limit, and figure 0 confirms that 
the minimum cross section approaches 2n at low couplings. 















Figure 3: Quantum absorption cross sections in the A c rs limit. The plot 
shows the absorption cross section as a function of energy, for small couplings, 
a = Mm/m p 2 <C 1. In the low-energy region plotted here, the wavelength of the 
Dirac particle is large compared to the black hole event horizon. Absorption is 
dominated by the lowest j = 1/2 angular momentum states (k = 1,-1 states), 
and the low-energy approximation of Unruh {2| is valid (see text). As a —■> 0, the 
minimum of the cross section approaches 2tt (dotted line). In the high energy 
limit, all cross sections return to the photon limit, 277 t. 
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6 Discussion 


We have shown that the absorption cross section for a Dirac wave in a classical 
Schwarzschild background can be calculated by matching ingoing solutions at 
the horizon to appropriate asymptotic forms at infinity. The analysis proved 
particularly simple in the Painleve-Gullstrand metric, though we stress that the 
cross sections calculated here do not depend on the particular choice of gauge. 
Furthermore, antiparticle solutions are be generated from particle solutions by 
the transformation 


(til, U 2 , E, k) i > (u\,u%, —E*,—k ). (33) 

It follows that the absorption cross section is invariant under charge conjugation. 

In the large-a limit (Rs A) we find that the cross section approaches the 
classical prediction of equation ©• When a ~ 1 (Rs ~ A) we observe energy- 
dependent oscillations about the classical value in the a-vs-E plot (figure Q. 
Oscillations of this nature were previously found by Sanchez for the massless 
scalar wave. The form of the oscillations depends on m, so represents a quantum- 
mechanical violation of the equivalence principle (figure Ell- 

In the low-energy regime, the j = 1/2 cross section 1311) given by Unruh 
is an excellent fit to our numerical results. In the a —> 0 (Rs -C A) limit we 
see that the minimum-possible cross section approaches 2tt(GMc~ 2 ) 2 . In the 
high-energy limit, all cross sections eventually converge on the photon limit of 
27 t t(GMc~ 2 ) 2 . 
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